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The Schrodinger equation for stationary states in a central potential is studied in an 
arbitrary number of spatial dimensions, say q. After transformation into an equivalent 
equation, where the coefficient of the first derivative vanishes, it is shown that in such 
equation the coefficient of r~ 2 is an even function of a parameter, say A, depending on 
a linear combination of q and of the angular momentum quantum number, say /. Thus, 
the case of complex values of A, which is useful in scattering theory, involves, in general, 
both a complex value of the parameter originally viewed as the spatial dimension and 
complex values of the angular momentum quantum number. The paper ends with a proof 
of the Levinson theorem in an arbitrary number of spatial dimensions, when the potential 
includes a non-local term which might be useful to understand the interaction between 
two nucleons. 
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1. INTRODUCTION 

Although the general framework of non-relativistic quantum mechanics is by now well 
known [1,2], while its predictions have been carefully tested against observations [3], it 
remains quite important to understand whether some familiar problems are a particular 
case of a more general scheme. For this purpose, we here consider the Schrodinger equation 
for stationary states in an arbitrary number of spatial dimensions, say q. This topic has 
always attracted interest [4-7], and we will see that a simple but deep result is found to 
hold independently of the particular value of q. Section 2 is devoted to a pedagogical 
formulation of the eigenvalue problem. Section 3 derives in detail the reduction of the 
Schrodinger equation to a simpler form. General results on scattering states are described 
in Sec. 4, and the Levinson theorem is proved in Sec. 5 in an arbitrary number of 
dimensions and in the presence of a non-local term in the potential. 

2. STATIONARY STATES IN q DIMENSIONS 

As is well known, a potential on lZ q which only depends on the distance from the 
origin, say r, is called spherically symmetric, or central. On denoting by A the Laplacian, 
we are interested in its action on the domain D of C°° functions with compact support 
away from the origin: D = Cq° (1Z q / {0}) . Each square-integrable function on 1Z q may be 
viewed as a function of r and of q — 1 variables £ on the (q — l)-sphere S' 9-1 . The solution 
by separation of variables of the Schrodinger equation for stationary states leads to the 
consideration of functions in D which are finite linear combinations of products f(r)g(£). 
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Such a set is dense in the Hilbert space of square-integrable functions on 1Z q , because [4] 

L 2 {U q ) = L 2 (TZ+, r^dr) <g> L 2 (S q -\ dQ), (2.1) 

where <ifi is the integration measure on S^ -1 . Equation (2.1) tells us that square-integrable 
functions on 1Z q belong to the tensor product of the space of square-integrable functions on 
the positive half-line 1Z+ (the measure being r q ~ 1 dr) with the space of square-integrable 
functions on S q ~ 1 . On functions of the form f(r)g(£), the operator —A acts by [4] 

(-A)/(r)<?(£) = (-^ - f(r)9(t) - ±f(r)B q _ ig (^ (2-2) 

where B q _i is the Laplace-Beltrami operator on L 2 (S q ~ 1 ) [8]. The operator is 
found to be negative-definite, with only point spectrum of finite multiplicity, and with C°° 
eigenfunctions. The Z-th eigenvalue, ki, is given by [4,8] 

h = -1(1 + q -2). (2.3) 

Thus, the Schrodinger equation for stationary states [5]: 



h 2 



A +U(r) *(f) = EV(r), (2.4) 



2m 

leads to the following equation for the radial part, ifi(r), °f the wave function [5,7]: 
d 2 (q-1) d 2m l(l + q-2)~ 



dr 2 r dr h r l 
on expressing the wave function as the product [5] 



V>(r) = 0, (2.5) 



$(r1^(r)#,..,Vi). (2.6) 
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The function g in Eq. (2.6) is known as the generalized spherical harmonic [5,8], and 
further indices have been omitted for simplicity of notation. Of course, the full solution 
is eventually obtained after summing all contributions of the kind (2.6), the summation 
being taken over the whole set of quantum numbers corresponding to the given value of q. 

3. REDUCTION TO A SIMPLER FORM 

One now wants to map Eq. (2.5) into an equation where the coefficient of the first 
derivative vanishes, so that a problem involving the Hilbert space of square-integrable 
functions on the positive half-line 7Z + with respect to the measure r q ~ 1 dr is mapped 
into a genuinely "one-dimensional" problem, in that one considers the square-integrable 
functions on the positive half-line 1Z+ with respect to the measure dr. For this purpose, 
one defines a new function y by means of the equation (see page 194 of Ref. [9]) 



where (3 is, for the time being, an unknown parameter. The request that, upon insertion 
of Eq. (3.1) into Eq. (2.5), one should get a second-order differential equation for y where 
the coefficient of ^ vanishes, identifies j3 as the solution of the algebraic equation 



tfj(r) = r y{r), 



(3.1) 



2/?+(<z-l) = 0, 



(3.2) 



and hence the desired transformation reads, eventually, 



i^{r) = r 



(g-i) 



y(r). 



(3.3) 
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On setting 



2mE 



(3.4) 



K 



2 ' 



h 



V(r) 



2m 



(3.5) 



ti- 



the resulting differential equation for y reads 



dr 2 



( 



/(/ + q - 2) + -{q 2 - Aq + 3) J ^ + k y{r) = V(r)y(r). 



(3.6) 



So far, what we have done is quite standard. However, unlike the treatment in Refs. [5,6], 
where the coefficient of r~ 2 is expressed in a form proportional to (p — l)(p — 3), with 
p = q + 2/, we would like to re-express it in a form which is an even function of another 
suitable parameter, say A, because this is indeed the case when q = 3, and one wonders 
what happens if the number of spatial dimensions is greater than 3. For this purpose, we 
remark that the familiar technique of completing the square leads to 




(3.7) 



and hence one finds 




(3.8) 



Thus, the desired parameter, A, should be defined as 



A = J+i(g-2), 



(3.9) 



5 



so that Eq. (3.6) may be cast in the very convenient form 



d 2 „ (A 2 -|) 



dr 2 



+ K- 



-V(r) 



y(r) = 0. 



(3.10) 



It is now clear that, upon performing the transformation (3.3), one obtains a second- 
order differential equation for y in complete formal analogy with what happens in TZ 3 . In 
particular, for scattering states, one has E > 0, and it is more convenient to define 



k 2 = 



2mE 



h 



2 ' 



(3.11) 



so that Eq. (3.10) reads 



— + k 2 

dr 2 



(A 2 -i) 



V(r) 



y{r) = 0. 



(3.12) 



4. GENERAL RESULTS ON SCATTERING STATES 

To appreciate the utility of Eq. (3.12), let us recall what one learns from the analysis 
of the partial- wave boundary conditions at r = and at infinity [10]. The point r = is 
a Fuchsian singularity [11] of Eq. (3.12), which admits two fundamental solutions, say ip\ 
and (f2i such that [10] 

^(r) = r A+ ^(l + o(l)) as r -> 0, (4.1) 
(p 2 {r) = r" A+ ^(l + o(l)) as r -> 0. (4.2) 

Since the operator in Eq. (3.12) is even in A, the fundamental solutions cpi and <p2 exchange 
their roles under the substitution A — > —A. It is hence convenient to write <p(\, k, r) and 
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ip(— A, k,r) in place of (pi(r) and <p2(r), respectively, where the parameter A is allowed to 
be freely specifiable in the complex plane [10]. 

Moreover, as r — > oo, one looks for the Jost solution /(A, fc,r) of Eq. (3.12), which 
satisfies the asymptotic condition [10] 



lim e ikr f(\,k,r) = 1. (4.3) 



A deep relation exists between the boundary conditions at r = and at r = oo [10]. In the 
former case, one finds that (p(X, k,r) is entire (i.e. analytic in the whole complex plane) in 
k 2 and analytic in the domain Re(A) > 0. If, for some positive C, 

\V(r)\ < Cr~ 2+e , (4.4) 

the analyticity domain can be extended up to Re(A) > —\e. The Wronskian W satisfies 
the property 

W v?(A, k, r), (p(-X, k, r)] = -2A, (4.5) 
and, if V is real-valued, the following Hermiticity condition holds: 

"^(A,fc,r)]* = V (A*,fc*,r). (4.6) 

In the latter case, one finds that f(X,k,r) is entire in A 2 and analytic in the domain 
Im(fc) < 0. If the potential satisfies the inequality 

\V(r)\ < Ce~ ur , (4.7) 
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the analyticity domain of / can be extended up to Im(/c) < ^. The Wronskian yields 



W f(X,k,r)J(X,-k,r) =2ik, 



(4.8) 



and, if V is real-valued, the following Hermiticity condition holds: 



* 



f{\\-k\r). 



(4.9) 



In our (/-dimensional problem, the form of Eq. (3.12) shows that all these properties 



freely in the complex plane, and not merely /. For example, from what we said above, the 
analyticity domain of <p(\, k, r) is expressed by the condition 



In other words, the complex angular momentum formalism is actually a particular case of 
a more general scheme, where both I and q are freely specifiable in the complex domain. 
Interestingly, one can even keep / real and still obtain complex values of A, provided that 
the parameter q is allowed to take complex values. Had one formulated from the beginning 
the eigenvalue problem in q dimensions, while bearing in mind the properties of scattering 
theory in 3 dimensions, the above results would have been clearer to the physics community. 

Yet another interesting effect of the factor \{q — 2) emerges in the analysis of potential 
scattering when the potential V in Eq. (3.5) is assumed to admit a Laurent expansion in the 



still hold. However, by virtue of the definition (3.9), it is I + | — 1 which is allowed to vary 




(4.10a) 



or, in equivalent form, 




(4.10.6) 
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annulus r e]0, oo[, so that Eqs. (2.5), or (3.6) or (3.12) acquire non-Fuchsian singularities 
both as r — > and as r — > oo [7,12]. One then finds solutions of Eq. (2.5) in the form 
[7,12] 



oo 



E 

n= — oo 

The fractional part of the parameter 7 can be found by solving a transcendental equation 
involving a Hill determinant which is an even periodic function of [7] 

l = l+\{q-2). (4.11) 

In other words, the polydromy parameter receives a contribution from q, and one may well 
allow q to take complex values at this stage. 

5. LEVINSON THEOREM WITH NON-LOCAL POTENTIALS 

The ground is now ready for studying a non-trivial application of the property ex- 
pressed by our Eq. (3.12). Following Ref. [13], we consider the stationary Schrodinger 
equation in a central potential when a non-local term contributes to such potential. This is 
motivated by what one knows about the interaction between two nucleons, which is purely 
local at sufficiently large distances, but acquires a non-local nature as the two nucleons 
approach each other. The novel feature in our analysis, with respect to Ref. [13], is the ar- 
bitrary number of spatial dimensions, so that, for scattering states, we study the equation 
(cf. (3.12)) 



d 2 , , 2 (A 2 - \) 



dr 2 



+ k — ± ^ - V(r) 



r 



y{r) 



r <?2 } J U(r,r')y(r')r' (q 2 } dr' . (5.1) 



In Eq. (5.1), the kernel U(r,r'), also called, with little abuse of notation, non-local po- 
tential, is taken to be real-valued, continuous, symmetric, and vanishing at large distances 
[13,14]: 

U{r,r') = U{r',r), (5.2a) 



U(r,r') ~ 0(r _1 ) as r -> 0, 



U(r,r') = Vr > r , 



(5.26) 
(5.2c) 



where r is a sufficiently large value of the radial coordinate (as we said before). The local 
part, V(r), of the potential, is taken to be real- valued, continuous, and such that 



V(r) ~ 0(r _1 ) as r -> 0, 



(5.3a) 



V(r) = Vr > r . 



It is also useful to introduce a real parameter, say for which one can write 



(5.36) 



V(r, fx) = (xV(r), U{r, r' , fx) = (J,U(r, r'). 



(5.4) 



The idea is that, as \i ranges from through 1, the rescaled potentials V(r, fj,) and U (r, r', fi) 
range from through the original values V(r) and U(r,r'). The radial equation (5.1) is 
hence replaced by 



dr 2 



+ k 2 - — — tt^- - V(r, jj) 



yk,\i r i v) 



J t/(r,r / ,^)y fc)A (r / ,^)r /(9 2 1) rfr / . 



(5.5) 
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We now consider Eq. (5.5) for two different values, say k (with solution yk,x) and k (with 
solution Vk x)-, and multiply the equations for yk,x(r, n) and y kX (r,ii) by y kX (r,y) and 
Vk,\{ r i AO > respectively. On taking the difference of the resulting equations, one finds 



d_ 

dr 



(yk,x-^y k ,x - y k ,x^yk,\^ +(k 2 - k 2 ^y k}X y K 



= r (q i * yk,X / U(r,r',v)y kX (r',n)r' (9 * ' ' dr' 



r 2 Vk,x U(r,r ,/J>)y k ,x{r ,y)r 2 dr. 



(5.6) 



Since, by regularity, both y kjX and y kX have vanishing limit as r 
(5.6) over the interval [0, ro] yields 



d_ d 



<9r' 



0, the integration of 



-2 \ Z" 7 " 



r (q 2 1) r > y kjX ( r ', n)y k Ar, n) U (/, r, /i) - U (r, r', ^) I dr dr' = 0, (5.7) 



where the right-hand side of (5.7) vanishes by virtue of the symmetry condition (5.2a). 
Since k ^ k by hypothesis, we can multiply both sides of (5.7) by -, A= — -, 7^2^ — , 



which yields 



(k 2 - k 2 ) 



1 d 

Vk,x( r o^) dr 



Vk,x( r > »)\ 



_ 1 d 

r ° y k , x (r , /j) dr 



Vk,x{r, fj)\ 



ro 



r ° y k ,\{r',ij) y k ,x(r',v) ^ 

yk,x( r o^)yk,x( r o^) 



(5.8) 
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At this stage, on taking the limit of both sides of (5.8) as k — > k, one finds 



d_ 

dE 



1 d 

VE,x(r, n) dr 



pro 

= -Ve 2 \^ y) I y 2 E,\( r ', y)dr' < o, 

r=r a JO 



(5.9) 



where the subscript k for y(r, /i) has been replaced by E = ^k , which is more convenient 
from now on. Similarly, one finds 



d_ 
~dE 



1 d 

VE,\(r, n) dr 



VE,\{r, (i) 



POO 

= y- 2 x (r , n) / y|, A (r r , rfdr' > 0. (5.10) 

+ Jr 



Equations (5.9) and (5.10) tell us that, as energy increases, the logarithmic derivative of 
the radial function at r$ decreases monotonically, whereas that at r$ increases monoton- 
ically. This is an expression of the Sturm-Liouville theorem in q dimensions for non-local 
potentials in a central field (cf. Ref. [13]). 

The matching condition at ro for the logarithmic derivative of the radial function is, 
of course, 



A X {E, fj) 



1 d 

Ve,\{t, n) dr 



VE,x(r, ft) 



1 d 



VE,\{r, (i) 



(5.11) 



The form (5.5) of the stationary Schrodinger equation is indeed quite involved, but some 
limiting cases are easily dealt with. For example, for a free particle, which corresponds to 
a vanishing value of \i, one finds ( J\ being the Bessel function of first kind and order A) 



VE,x(r, 0) = \l -nkr J\(kr), 



(5.12a) 
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when E > 0, and 



(5.126) 



with «= vEpl 5 if £■ < 0. 

In the interval [ro, oo[, both 1/ and t/ vanish, and for positive values of E two oscillating 
solutions of (5.5) exist, so that the general solution reads (cf. Ref. [13]) 



VE,\{r, fi) 



-nkr 



J\(kr) cos T] X {k,n) - N x (kr) smrj X (k, fi) 



(5.13) 



where r]\(k, fx) is the phase shift, and N\ is the Neumann function of order A. The matching 
condition (5.11) leads to a very useful formula for the phase shift, upon using the result 
(5.13), i.e. 



tanr]x(k,n) = 



J\(kr ) 


A x (E,/j) 


]„J' x (kr ) 
J\{kr ) 


1 

2r _ 


Nx(kro) 


Ax(E,fi) 


h N' x (kr ) 


1 

2r _ 



(5.14) 



Equation (5.14) provides the key tool for proving the Levinson theorem, jointly with a 
careful analysis of matching conditions. Such a theorem leads to a relation between the 
total number of bound states and phase shifts at zero momentum. Here, following Ref. 
[13], we shall agree that the phase shift is determined with respect to the phase shift 
r)\(k,0) for a free particle, where, by definition, one chooses r) X (k,0) = 0. With such a 
convention, the phase shift is determined completely as \i increases from to 1. 
If E < 0, the only square-integrable solution of (5.5) is 



Ve,x 



D i(X + l)7T/2 



ttkt H^' (inr) 



(5.15) 
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where is the standard notation for Hankel functions of first kind and order A. Thus, 
the right-hand side of the matching condition (5.11) reads 



d 



yE,\(r,Q) dr 



iKH^\mro) 1 
n= r +~ fff^ro) 2r ' 



(5.16) 



As E — > , the right-hand side of (5.16) reduces to 



-A+- ) — =px, 
2/ r 



(5.17) 



whereas it tends to — k as E — > — oo. Moreover, the solution (5.12b) satisfies the condition 
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2/E,A(r,0) <9r 



ys,A(r, 0) 



mJ^(mro) 1 



(5.18) 



whose right-hand side tends to 



A + 7: ) — = pa as E 
2/ r 



0" 



(5.19) 



Thus, no bound state exists when fx = 0. Unlike the two-dimensional case studied in Ref. 
[13], not even the so-called "half-bound states" (for which p\ = p\) may occur, because A 
can never vanish if / is > and q is bigger than 2 (see (3.9)). 

Another crucial remark is that, if A\(0, p) decreases across the value p\ as p, increases, 
an overlap between two ranges of variation of the logarithmic derivative on the two sides 
of ro occurs. Bearing in mind the Sturm-Liouville theorem expressed by (5.9) and (5.10), 
such an overlap means that the matching condition (5.11) can only be satisfied by one 
particular value of the energy, and hence a scattering state is turned into a bound state. 
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In general, as p increases from to 1, each time A x decreases across p x , a scattering state is 
turned into a bound state for the above reasons. By contrast, each time A\(0, p) increases 
across p\, a bound state is turned into a scattering state. The number of bound states 
is then equal to the number of times that A x (0) decreases across p x as p ranges from 
through 1, minus the number of times that A\(0) increases across the value p\. The next 
task is now to prove that this difference equals r)\(0), the phase shift at zero momentum, 
divided by n. 

For this purpose, we have to evaluate tan p) when k « because 

t/a(0,/i) = \im m(k,ii). 

By virtue of the exact result (5.14) and of the limiting behaviour of Bessel functions at 
small argument, one finds, to lowest order in kro, 

taxir}x(k, p) 



n(kr ) 2X 


A x (0,p)-(\ + 


2>r 


2 2A Ar 2 (A) 


A x (0,p) + (\- 


-)- 

1) r _ 



7i (kro 



,2A 



A x (0,n) + (j> x -±) 



22A A r 2 (A) [A x (0,p)-p x ] 



(5.20) 



where we have neglected, for simplicity, higher-order terms in kro i n the denominator 
(strictly, higher-order terms should be included in the small-(/cro) expansions, whenever 
one has exact cancellation of leading terms). Unlike Eq. (21) of Ref. [13], our Eq. (5.20) 
expresses the tangent of the phase shift in terms of a formula involving the T-function, 
because A is not necessarily an integer. However, it remains true that tanr] X (k, p) tends 
to as k — > (since A > 0), and hence r)\(0, p) is always equal to an integer multiple of 
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7r. This means that the phase shift changes discontinuously. Besides, the exact formula 
(5.14) can be used to prove that the phase shift increases monotonically as the logarithmic 
derivative A\(E, y) decreases [13]. Thus, rj\(0,y) jumps by n if, for k sufficiently small, 
tanr]\(k, y) changes sign as A\(E, y) decreases. To sum up, when y ranges from through 
1, whenever A\(0, y) decreases from near and larger than the value p\ to smaller than p\, 
the denominator in (5.20) changes sign from positive to negative, leading to a jump of 
rj\(0, y) equal to it. By contrast, whenever A\(0, y) increases across p\, rj\(0, y) jumps by 
—tv. Bearing in mind what we said in the paragraph after Eq. (5.19), we conclude that 
77a (0) divided by ix is indeed equal to the number n\ of bound states [13]: 

Vx (0)=n x n, (5.21) 

which is the form of the Levinson theorem with an arbitrary number of spatial dimensions 
in the presence of non-local potentials (see (5.2o)-(5.2c)). 

Needless to say, the analysis in the present section relies heavily on the work in Ref. 
[13], but provides a non-trivial application of Eq. (3.12) and also a non-trivial generaliza- 
tion of the work in Ref. [13]. There is therefore increasing evidence that the formulation 
of quantum-mechanical problems in an arbitrary number of spatial dimensions may lead 
to new developments in scattering theory with complex or real values of the parameter 
A defined in Eq. (3.9) (whereas, in the past, more emphasis had been put on large-g 
expansions in quantum mechanics [5]). 
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